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Abstract 



The quantum clock of Salecker-Wigner-Peres is used, by performing a post-selection of the final state, to obtain average 
transmission and reflection times associated to the scattering of localized wave packets by static potentials in one 
I dimension. The behavior of these average times is studied for a gaussian wave packet, centered around a tunneling wave 
number, incident on a rectangular barrier and, in particular, on a double delta barrier potential. The regime of opaque 
£N) barriers is investigated and the results show that the average transmission time does not saturate, showing no evidence 
£j-jOf the Hartman effect (or its generalized version). 

Keywords: Salecker-Wigner-Peres quantum clock, time dependent tunneling, average tunneling times, Hartman effect 



< 

in 



43 
i 

■4— > 

c 



> 

m 
O 

od 
O 



> 



1. Introduction 

As it is well known, the Hamiltonian of a quantum 
mechanical system must be bounded from below due to 
stability requirements and, as a consequence, there is no 
self-adjoint time operator canonically conjugated to the 
Hamiltonian p] (see, however, [2]), hence posing serious 
difficulties to define a meaningful time scale for quantum 
phenomena such as tunneling. 

Several approaches have been considered to circumvent 
the above difficulty, leading to many definitions of quan- 
tum tunneling time, among which the dwell time and the 
phase time are the most prominent (see [3J |H [3] and 
references therein). However, the time scales proposed un- 
til now are generally plagued by difficulties which prevent 
their interpretation as a traversal time in the general case. 
For instance, a serious controversy in the interpretation of 
many such time scales is related to the Hartman effect, 
which asserts that in the stationary case the time satu- 
rates for a static opaque potential [7] (also see [El IS] for a 
generalized version). 

However, one cannot properly address the problem of 
quantum tunneling time in the stationary case because 
the very concept of a particle having a definite tunneling 
time requires that it be represented by a localized wave 
packet (see, for example, [TO])- With this in mind, in this 
work we reexamine the tunneling of a wave packet through 
a static potential by using the quantum clock formalism 
introduced by Salecker and Wigner [TT] and further re- 
fined by Peres [1 2 . Our main concern here is to scrutinize 
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whether the Hartman effect, both in its original and gen- 
eralized versions, emerges or not for a wave packet. 

The Salecker-Wigner-Peres (SWP) clock QTJ [12] (see 
also [T3J1H]) consists of a quantum rotor weakly coupled 
to the system (the tunneling particle) , and it can be shown 
to be equivalent to the Larmor clock [TU] US US] , prepared 
and observed in a particular set of states, in the limit of 
large spin - for a detailed treatment see [17] . In this frame- 
work a "time operator" (not canonically conjugated to the 
Hamiltonian) can be defined |12j . and the application of 
such an operator to quantum tunneling was extensively an- 
alyzed by Leavens and McKinnon [18], who found that its 
expectation value only produces a meaningful result after 
a calibration procedure is adopted for the clock readings. 
Recently, Park [E5] considered the one-dimensional scat- 
tering of a wave packet by a rectangular potential barrier 
and found that the expectation value of that "time op- 
erator" in the asymptotic final states gives the dwell time 
averaged with respect to the components of the initial wave 
packet, which, however, cannot be considered a traversal 
time because it does not distinguish between transmitted 
and reflected particles [19] , 

In this paper, instead of associating the clock readings 
to expectation values of such a "time operator", we fol- 
low Peres' approach in the study of the time of flight 
[T2"] . According to this method the transmission (reflec- 
tion) time, which is obtained by following the peak of the 
clock's wave function, is given by the time derivative of 
the phase shift of the transmission (reflection) amplitude 
with respect to the perturbation potential. We argue that 
in order to address the question of the time spent in the 
barrier region by a particle that is eventually transmitted 
(reflected) one must perform a post- selection of the final 
state. Then, making use of an asymptotic condition ensur- 
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ing that for asymptotic times the incident and the trans- 
mitted/reflected wave packets are well localized far from 
the potential, it is possible to trace out the particle's de- 
grees of freedom in order to "read" the clock's final state. 
The average times obtained through this procedure are 
shown to have sensible physical properties in the context 
of quantum tunneling of gaussian wave packets through 
static symmetric potentials. 

This work is organized as follows. In section[2]we briefly 
review the formalism introduced by Peres [T2] and its 
generalization to the time dependent situation (also see 
[H1[H]) > introducing the necessary steps to take into ac- 
count only particles that are eventually transmitted (re- 
flected) by the potential. We discuss some properties of 
the average times obtained and compare them with other 
time scales suggested in the literature. In Section [3] we 
study the behavior of these average times for the scatter- 
ing of a gaussian wave packet, centered around a tunneling 
wave number, by a single rectangular barrier and a double 
delta barrier, with emphasis on the opaque regime and the 
possible emergence of the Hartman effect and, in partic- 
ular, its generalized version. In section [4] we present our 
conclusions. 



2. The Salecker-Wigner-Peres Clock 

The Salecker-Wigner clock jTTj , as considered by Peres 
|12) . is just a quantum rotor with Hamiltonian given by 
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with to = 2w/(Nt), where r is the time resolution of the 
clock and N = 2j + 1 is the dimension of the Hilbert space 
(with j being a non-negative integer or half- integer) . 

The eigenfunctions u m {9) = ^=e" n0 (0 < 9 < 2vr) of 
H c have eigenvalues V m — mhaj with m = —j, . . . , +j. A 
more convenient orthogonal basis for the clock is given by 
v r (6) (r = 0,...,N-l) m 



V r (9) 



m— ] u m \a) , 



N 



(2) 



which has the advantage of having a sharp peak at 9 = 
2irr/N for large N, so that the clock's pointer indicates 
the rth "hour" (with uncertainty ±^)- Under time evolu- 
tion such states change as v r (9) — > v r (9~ oot) , thus making 
possible to associate a time duration t to an arbitrary in- 
terval just by observing the translation of the peak of a 
free clock's wave function v r in that interval. 

We are interested in measuring the time associated with 
the one-dimensional scattering of a particle of mass /i by a 
(real) static potential V(z) confined to the region < z < 
a. Following [12], the particle and the clock are coupled 
through the interaction Hamiltonian 



H I =T{z)H c , 



(3) 



where the projection operator V(z) satisfies V(z) = 1 if 
z € [0, a] and equals zero otherwise. The clock is assumed 
to be initially at the state vq(9). Thus, the general form 
of the initial wave function of the whole system is 



\T/. 



1,z,t) = $(z,t)v (9), 



(4) 



where $(z,i) is a wave packet corresponding to a parti- 
cle incident from the left of the potential which, from the 
symmetries of the problem [2D], can be written as 



dk 



2?r 



(5) 



with all integrals in k going from — oo to +oo. The Fourier 
coefficient A{k) is assumed to be centered around a posi- 
tive k$ corresponding to a tunneling energy. 

After the scattering the particle and the clock become 
entangled due to the interaction Hi and the general form 
of the wave function is [13] 



where 



dk 
2^ 



A(k)ip {m) (k,z)e- 



(6) 



(7) 



The wave function ^>(9,z : t) satisfies the Schrodinger 

2 

equation with the Hamiltonian H = |- + V(z) + Hi. The 
orthogonality of the u m (6)'s, together with the fact that 
they are eigenfunctions of H c , imply that ip^ m ^{k,z) sat- 
isfy the time-independent Schrodinger equation for energy 
E(k) and potential 



V^ m \z) = V{z)+V m V{z) 



(8) 



where V m = mhuj. In order to ensure a weak coupling 
between the clock and the system it is assumed that V m is 
a small perturbation, i.e., \E\ 3> V m and |V(,s) — E\ ^> V m . 

Since we are interested in measuring the transmission 
time asymptotically after the interaction, let us focus on 
the transmitted wave packet. The solution of the above 
stationary problem in the region to the right of the poten- 
tial V^ n \z) is 



'tPtransik, z ) 



T {m \k)e ikz 



i^T (k)+ikz 



By using the weak coupling condition we can assume that 
\T( m )(k)\ ~ \T(k)\ and expand the phase <p${k) up to first 
order in V m H2l: 



v ^ { k)^ VT {k)+V m (-^-^\k) 



(9) 



Vm=0 



where (firik) denotes the phase of the transmission coef- 
ficient corresponding to the unperturbed potential V(z). 
Finally, substituting the above expressions in (|6|-([7]), and 



using the explicit form of the functions u m (9), the asymp- 
totic transmitted wave function is given by 



dk 
2^ 



A(k)T(k)e^ kz -^ 
xvo(6-u%(k)) , (10) 



where 



(ii) 



is the transmission time used by Peres in [T2]. A similar 
expression can be obtained for the reflection time tf-(k) 
in terms of the reflection phase Lp^\k). It can be proved 
that these clock times satisfy [T3] (also see [2"T] ) 



r D (k) = \T(k)\ 2 ti (k) + \R(k)\ 2 tf(k), 



(12) 



where ro(fc) stands for the stationary dwell time 15], and 
T(k) and R{k) are the transmission and reflection ampli- 
tudes, respectively. It is important to stress that, for the 
clock times used here, the above relation can be proved to 
be a direct consequence of Schrodinger's equation. Thus 
the criticism to such a relation for supposedly not being 
compatible with quantum mechanics due to the absence of 
an interference term, which sometimes appear in the lit- 
erature (see, e.g., [4]), does not apply to the present case 
(for details see Q3]). 

It should be noticed that, in order to strictly enforce the 
conditions \E\ ^S> V m and \V(z) — E\ ^> V m , the integral 



in ( 10 ) should in principle be truncated by a small cut- 



off excluding the regions in which these conditions are not 
warranted. However, the continuity of the integrand allow 
us to extend the integration over the whole fc-space up to 
first order for a very small V m , consistent with the approx- 
imations in @. This requires a large time resolution r for 
the clock, since it is related to the reciprocal of V m [13] . 
However, as observed by Davies [22]; there is no difficulty 
since the clock can be treated within the weak measure- 
ment theory of Aharonov et al. [23] [21] and a precise value 
can be obtained for a large enough ensemble. 

Here we are interested in finding an average transmission 
time. Therefore, we must consider only the sub-ensemble 
of all scattered particles that are eventually transmitted, 
i.e., we must perform a post- selection of the final state (in 
addition to the preselection of the initial one). The post- 
selected (transmitted) state can be described, with a slight 
change in notation, by the density matrix of the particle- 
clock asymptotic state g trans = a trans) (^transl, where 
a is a normalization constant and |$ trans) corresponds to 
the state (10). In order to observe the final state of the 



clock we assume an asymptotic condition, stating that long 
after the interaction the transmitted wave packet is com- 
pletely localized to the right of the potential. Then, we 
can trace out the particle's degrees of freedom (by extend- 
ing the z-integral over the whole space) and obtaining the 



clock's reduced density matrix for the transmitted sub- 
ensemble, given by 



trans 



^\A(k)T(k)\ 2 
x \vo (0 - ut T c {k))) {v Q (e - ut T c {k)) | , (13) 



trans \ 



= 1 im- 



where the normalization requirement Tr(p 

plies a={J§\A(k)T(k)\ 2 y 1 . 

Now, since the clock's reading is associated with the 
peak of its wave function, we define the average transmis- 
sion time as (t£) — Tr (^g^ Tans ^j , where 9 is an operator 
giving the peak of Vk [2S]- Then, 

(O = ( dkp T (k)t*(k), (14) 



is the probability density to 



i n\ \A(k)T(k)\ 2 

where p T (k) = f ^l4k)T(k)\» 
find mode k in the final transmitted wave packet. An 
analogous procedure for the reflected sub-ensemble results 
in the average reflection time: 



with pn(k) ■ 



dk PR (k)t^(k) 



\A(k)R(k)\ 2 



(15) 



fdk\A(k)R(kr the corres P ondin g Probability 



density to find mode k in the reflected wave packet. 

In the special case of a very sharply peaked wave packet, 
|A(fc)| 2 /(27r) ~ 8(k - k ), definitions (pl-iTEl reduce to 



the corresponding stationary times, as expected. However, 
for more physical incident wave packets, which necessarily 
involve a finite dispersion in wave numbers, these defini- 
tions have the important property of weighting the contri- 
bution of each tj(k) (t^(k)) by the probability density of 
having the mode k in the final transmitted (reflected) wave 
packet, thus minimizing the contributions of modes leading 
to the most improbable results. This is in contradistinction 
with the approach of considering the density matrix of the 
transmitted and reflected wave packets together, which is 
equivalent to averaging the dwell time with respect to the 
initial wave packet components - a procedure that leads 
to many well known difficulties [HJ [21] . 

The expressions ( 14 )-( 15 ) coincide with the real part 
of the complex time obtained by Leavens and Aers 16J 
in their analysis of double barrier resonant times using 
the Larmor clock (see also [2Z] which deals with similar 
averages for the time obtained in [2E])- These definitions 
could also be inferred from the analysis of the Larmor clock 
by Falck and Hauge [10] . The main advantage of the above 
treatment (in which only the position of the clock's pointer 
is observed [211 122]) is that it introduces a unique real time 
scale, thus avoiding the problem of interpreting complex 
times. 

The above average reflection and transmission times can 
be used to generalize relation ( 12 ) for a wave packet. First, 



we observe that the asymptotic probability to find the par- 
ticle in the region at the right of the potential (i.e., the to- 
tal probability for the particle to be transmitted) is given 



3 



by P T = limt^oo^ 00 dz J Q * dd\% rana (0, z, t)\ 2 . Under the 
asymptotic condition stated above the z-integration can be 
extended to the whole space and this probability reduces 
to 



dk 

/V / -\A{k)T(kt 



(16) 



Analogously, assuming a similar asymptotic condition 
stating that long after the interaction the reflected wave 
packet is completely localized to the left of the poten- 
tial, the total probability for the particle to be reflected 
is P R = J ^\A(k)R(k)\ 2 . From the normalization of the 
initial wave packet these probabilities satisfy Pt + Pr = 1 
(as required, since long after the interaction transmission 
and reflection are complementary outcomes). Then, defin- 
ing the mean dwell time as [3^1 EH] 



TD 



dk 

2^ 



\A(k)\ 2 r D (k) . 



(17) 



where the overbar indicates averages over the initial wave 
packet, the above definitions, combined with (12 1, imply 
that 



td 



1 



dk 
2^ 



\A(k)T(k)\X (k) 



dk 
27 



\A(k)R(k)\ 2 t?(k) 



P T {t T c )+P R {t R ) 



(18) 



This relation says that the usually defined average of the 
dwell time over the initial wave packet (which does not 
take into account whether the particle is eventually re- 
flected or transmitted) can be viewed as an average over 
the transmission and reflection average times, weighted 
by the total probabilities for transmission and reflection, 
respectively. It must be stressed that a necessary require- 
ment for this interpretation is that the asymptotic condi- 



tions be satisfied. Relation (181 appears as the real part of 
a similar complex relation in |16j - again, the advantage 
here is that the SWP clock provides only real times. 



3. Tunneling through symmetric barriers 

In this section we study the average transmission time 
for a gaussian wave packet incident from the left on two 
kinds of symmetric localized potentials, namely a single 
rectangular barrier and the double delta barrier. In both 
cases, we shall be mainly concerned with the behavior of 
(tl) in the regime of opaque barriers (characterized by very 
small transmission amplitudes for non resonant tunneling 
energies), in order to examine the emergence or not of 
the Hartman effect [7] (single barrier) and its generalized 
version [5] (double delta barrier). 

We assume that at t = the particle is in a state de- 
scribed by a right-moving gaussian wave packet centered 



around a tunneling wave number k > 0, and spatially 
centered around a position z < 0, i.e., 



$(z,0) 



1 



(2tt)V4, 



exp < ikoz - 



{z-z f 
Aa 2 



(19) 



where a is a scale measuring the spatial extent of the ini- 
tial wave packet (its inverse is proportional to the wave 
packet's momentum dispersion). 

In what follows we will always make use of the fact that, 
as a consequence of (12), for symmetric potentials t£(k) 



and tj(k) coincide with the dwell time To(k) [10]. There- 



fore, the average times ( 14 ) and ( 15 ) are just the dwell time 
averaged over the transmitted and reflected sub-ensembles, 
respectively. 

3.1. The single rectangular barrier 

Let us now consider a square potential barrier of height 
Vq and width a, V(z) — VoQ(z)Q(a — z). For a tunneling 
incident energy, E < Vq, the stationary dwell time is well 
known and given by |15) 



r D (k) = 



2fi k 



[(k 2 + q 2 ) tanh(ga) + qa(q 2 - fc 2 )scch 2 (<7a)] 
[Aq 2 k 2 + (q 2 - k 2 ) 2 tanh 2 (ga)] 



(20) 



where k = j-y/2fiE and q = jr^2fi(Vo — E). For energies 
above the barrier, E > Vq, the stationary dwell time can be 
obtained from the above expression through the substitu- 
tion q — > iki, with k\ = \^j2[i(E — Vq)), while the trans- 
mission amplitude T(k) can be found in any textbook on 
quantum mechanics (see, e.g., [20). In the opaque limit, 
characterized by qa ^> 1 for a given E < Vq, the station- 



ary dwell time ( 20 1 becomes independent of the barrier 



width a, which is the Hartman effect [7J |5] ■ This effect has 
been at the origin of an intense debate in the recent lit- 
erature concerning, among other problems, the apparent 
superluminality of the tunneling process (see, for exam- 
ple, [51 H [IB IH and references therein). It should be 
noticed that a time can be defined, based on an analogy 
with fluid mechanics |31j . such that in the stationary case 
it reduces to the dwell time divided by the transmission co- 
efficient and, therefore, does not saturate with the barrier 
width [3T] (also see [351 [33J [33] ) ; but a complete analysis 
of time dependent tunneling using such time scale is yet to 
be achieved [3X1135] . 

In Figure 1 we investigate the typical behavior of the av- 



erage times (14) and (15) with increasing barrier widths. 



For comparison we have also plotted the average reflection 
time ( 15 ), the mean dwell time (|17|), and the time spent by 



the incident packet (characterized by its peak wavelength 
fco) in the region [0,a] in the absence of the potential, i.e., 
tfree = ■ The choice of parameters in these figures 
ensures that the initial wave packet (19) is, to an excel- 



lent approximation, confined to the left of the barrier (at 
t = the probability to find the particle in the region 
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Figure 1: (Color online) Behavior of the average times with respect 
to the barrier width a. All quantities are expressed in atomic units 
(a.u.): jU = 1, a = 10, fcrj = 0.7, zq = —8<j, and Vq = 0.5. In- 
set: zoom of the region of small and moderate barrier widths [the 
horizontal dotted line indicates the saturated value of the stationary 
dwell time to(/co)]. The free time, tf ree , is plotted for comparison. 



z > is less than 10 -15 ). Also, tunneling components are 
dominant in the Fourier decomposition of the initial wave 
packet. We observe that in the extreme opaque regime (in 
this case characterized by very large barrier widths), while 
td saturates to a value independent of a, \t£) increases 
linearly with the barrier width. In the inset, showing in 
detail the region of small to moderate barrier widths, we 
observe that for thin barriers all the average times are 
very similar. The similarity between (iJT) and (tf') [hence, 
from (18), their similarity with fn] is expected, since for 
thin barriers the post-selected probability density distri- 
butions pa and pt are almost indistinguishable from the 
incident one, Pi nc . As the barrier thickness increases (t^) 
becomes smaller than the free time and soon saturates, as 
expected for a wave packet centered in a tunneling wave 
number, since one would anticipate that the wave function 
cannot penetrate beyond a certain barrier depth. The av- 
erage transmission time, (t^), also becomes smaller than 
the free time for intermediate values of the barrier width, 
indicating that the barrier speeds up the wave packet in 
the transmission channel; however, in the very opaque re- 
gion its behavior suffers an abrupt change, growing sud- 
denly until it attains a regime of linear growth - as shown 
for large a's in Fig. 1, in this regime (t£\ > tf ree (The 
wave packet width exerts an important role in determining 
the values of a for which this abrupt change occurs, with 
the barrier width for that change increasing with a; for a 
similar conclusion using a different time average see [36] ) . 



Therefore, we observe that (ij) never saturates, as one 
would expect of a well defined average transmission time. 
Finally, as it is well known, the mean dwell time initially 
increases with a and soon saturates (the Hartman effect); 
this is a consequence of the fact that in the opaque region 
Pr ~ 1 and the reflection term dominates the r.h.s of Eq. 



It is instructive to compare the qualitative behavior of 
the average times shown in Fig. [T] with the results pre- 
sented in [33], even though [33] considers the time inde- 
pendent situation. In both cases the reflection times sat- 
urate and the transmission times always increase with the 
barrier width, not showing any sign of the Hartman ef- 
fect. However, there are important differences: in 34J 
the reflection time starts from infinite for a — and it 
is very large for very small barrier widths before staurat- 
ing; furthermore, in |34j the transmission time tends to 
increase exponentially with a in the extreme opaque limit 
(also see [3U [321 1331 131] )■ These are unexpected properties 
and there is some debate about the interpretation of these 
times [3UI371H] (nevertheless, the transmission dwell time 
used in [34] has been shown to have practical use in the 
study of the half life of radioactive nuclei decaying by the 
emission of an a-particle [3H])- By contrast, (t^) is small 
for small barrier widths, a plausible result since this time 
is interpreted as the average time spent into the potential 
region by the wave packet which is eventually reflected. In 
addition, (t£) tends to increase linearly with a in the ex- 
treme opaque limit. Such good properties of (i[T) and {tf) 
emerge as a consequence of averaging over pos-selected 
sub-ensembles. 

The behavior of the average transmission time observed 
in Fig. 1 is a consequence of the fact that for the trans- 
mitted wave packet in the opaque regime (i<f) is domi- 
nated by above the barrier components. If one truncates 
wave packet to include only under the barrier components, 
the corresponding "truncated" transmitted wave packet 
is dominated by the contribution of a single component, 
completely lacking localizability in the opaque limit and 
approaching a plane wave and, therefore, saturating [but 



(181. 



in this case, strictly speaking, the definition ( 14 ) does not 
apply]. The above the barrier dominance is, of course, well 
known (by the way, methods to avoid such dominance, by 
compensating for the reflection of lower k components in 
the state preparation, have been developed and result in 
average times that saturate [55]). 



3.2. The double delta barrier 

The generalized Hartman effect states that for two po- 
tential barriers the tunneling time becomes independent 
of the distance between the barriers in the opaque limit 
[3 El HOI E] (see, however, [32] for an alternative view 
of the stationary case using a multiple peak decomposi- 
tion). In order to investigate this effect we will consider 
the double delta barrier V(z) = r yS(z) + r )5(z — d), where 
7 > gives the barrier strength and d is the separation be- 
tween the barriers and we will assume that the SWP clock 
runs only when the particle is in the region < z < d. 
This potential has the advantages of being the easiest way 
to study the generalized Hartman effect [33J and that all 
the wave packet components are tunneling (for a study of 
stationary traversal times for this potential, see [33])- 



5 



For a stationary wave incident from the left, the trans- 
mission coefficient is given by (in the limit V m = 0) 



d=10 



d=25 



\T{k)Y 



1 



1 + 4a 2 [asm(kd) + cos(kd)Y 



(21) 



jpr. In the opaque limit, characterized by 



where a 

a»l (equivalently, 7 k), |T(fc)| 2 goes to zero almost 
everywhere (as 7 -4 ), except at the (countable) infinite set 
of resonances k n . The resonances are given, for any 7 > 
and d > 0, by the solutions of the transcendental equation 



tan 



h*k 



(22) 



with the location of the nth resonance strongly dependent 
on d, in such way that the spacing between successive res- 
onances decreases with d (see [46] and references therein). 
For 7 — > 00 the resonance spectrum tends to the spectrum 
of an infinite rectangular well of width d, i.e. k n 
(withn = ±l,±2,±3,---) |2JJ]. 

The stationary dwell time in the region between the bar 



,1 



rieis, as follows from (11), is given by 
T D (k) 



t c (k) 



hk 2 



(l + 2a 2 )kd + 2a sin 2 (kd) - a 2 sin(2kd) 

X 2 * 

1 + 4a 2 [asin(fcrf) + cos(fcd)] 
and at the nth resonance, k n , its value is 



TD{k n ) = 



hkl 



,M7 

'h 2 



9 2 

2 M 7 +fc 2 



(23) 



(24) 



as can be found by substituting ( 22 ) into ( 23 ) . For a fixed 



d, the opaque limit of the resonant dwell time behaves 
asymptotically as To(fc n ) ~ ^s-^r d, thus diverging in the 
extreme opaque limit 7 — > 00. On the other hand, the 
off resonance dwell time behaves in the opaque limit as 
Tjy(k) cx 7 -2 , thus characterizing the generalized Hartman 
effect by the fact that Tjj(k) — > as 7 — > 00 (however, such 
a convergence to zero is non-uniform in any neighborhood 
of a resonance) . 



Let us now consider a gaussian wave packet sharply 



peaked around ko, eq. (19 1, and investigate the behavior of 
the average transmission time in the opaque limit (which is 
now characterized by 7 3> ko). It is useful to first analyze 
the wave-number distribution of the transmitted and inci- 
dent wave packets, given by the probability densities pr{k) 
andp mc (fc) = |A(A;)| 2 /(27r), respectively. As shown in Fig- 
ure [2j for small d there are few, if any, resonances in the 
domain of Pinc(k). However, as d increases the resonances 
become densely distributed within that domain and we 
expect the transmitted packet to be dominated by them. 
Such a dominance is especially pronounced in the opaque 
regime, where the off resonance components are strongly 
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Figure 2: (Color online) Wave number distribution for the trans- 
mitted (solid curve, non normalized) and the incident wave packets 
(dashed curve), for opaque barriers. Atomic units are used, with 
7 = 16, a = 6 and ko = 1.2. 
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Figure 3: (Color online) Behavior of the total probability for the 
particle transmission with increasing strength 7. We observe that 
this probability goes to zero as the barrier strength 7 increases, irre- 
spective the eventually large number of resonant components in the 
initial packet wave number distribution. In this plot ko = 1.2 and 
a = 6 in atomic units. 



suppressed (a similar behavior was found for the double 
square barriers in [46] ) . It is important to emphasize that 
even when the transmitted spectrum is dominated by res- 
onances, in the extreme opaque limit the wave packet total 
transmission probability tends to vanish, since in the limit 



7 — y 00 the integrand in ( 16 ) tends to zero almost every- 
where, except in a zero-measure set (the countable set of 
resonances) - this is illustrated in Figure [3] 

Figure^ displays (if), (t^) and td, in a logarithmic 
scale, for small to moderate barrier separations, and for 
7 in the opaque domain. One observes that the average 
reflection time only deviates significantly from the mean 
dwell time when resonances are present around fco, since 
there the reflection coefficient |i?| 2 is very small. Thus, 
when resonances are densely distributed around fco (e.g., 
when d is large) one expects significative differences be- 
tween (t^ and td, as it is shown in Fig 4p, for vary large 
d's. 

In Figure Hb we observe the typical behavior of the aver- 



age times for a wide range of values of the barrier separa- 
tion d. Contrary to the generalized Hartman effect state- 
ment, this figure reveals that (t;T) always varies with d 
(the plot corresponds to an initial wave packet having a 
narrow width Ak ~ - around ko) and, in fact, for large 
enough barrier spacing (ff) tends to grow linearly with d 
(although not shown here, we observed a similar behavior 
for several sets of values of 7 and a) . This is a consequence 
of the fact that for large d any small neighborhood of a 
given k will be densely populated by resonances (see Fig. 
[2]) and, as d varies, these resonances (of different orders) 
enter and exit that neighborhood, with the net effect of a 
linear dependence on d according to eq. ( |24| ) (with n not 
fixed, in a way that k n is always within such a neighbor- 
hood). On the other hand, the behavior of (i^T) for small 
and moderate barrier separations d tends to be oscillatory 
because there are few resonances in the transmitted spec- 
trum (see Fig. [3]) and, as d increases and a new single 
resonance eventually enters the domain of the transmitted 
spectrum, the peak of To{k) at that resonance causes a 
peak in (t^). 

Summarizing, Figures [4^, and |4Jd show that the aver- 
age transmission time depends on d, showing no evidence 
of the generalized Hartman effect. This is in opposition 
to the conclusions of 27] (using a different definition 
of average time, but with similar parameters) who assert 
that in regions of small d located between the peaks associ- 
ated to successive resonances (such that k n — k n -% 3> Afc) 
the generalized Hartman effect occurs, as characterized by 
plateaus in between the peaks in the plot time versus d 
(here such plateaus would be at zero) - no such plateaus 
appear for (t£) (in fact, not even for td) in the corre- 
sponding regions in Fig. [4Jl In addition, Figure |4]d shows 
that in the opaque regime 7 >> 1 the mean dwell time 
is dominated by the average reflection time, as one would 



expect from relation (18 1, and for very large d the mean 



dwell time tends to behave as the free time. In that range 
of d, (i^T) shows that the particle slows down for large 7, 
as expected due to the multiple reflections occurring in the 
region between the barriers. 



4. Concluding Remarks 

We considered the extension of the SWP clock to treat 
the one-dimensional tunneling of a particle, represented by 
a wave packet, through a static potential and addressed 
the question of how much time a particle that is even- 
tually transmitted (reflected) spends in the potential re- 
gion. Taking into account that to properly analyze this 
question one must consider only the sub-ensemble of the 
scattered particles that are finally transmitted (reflected) , 
which corresponds to a post-selection of the final state, 
the SWP clock attributes a well-defined average time for 
these process, namely eqs. (14 1 and (15). These are to 
be compared with the mean dwell time obtained in |19j 
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Figure 4: (Color online) Behavior of the average transmission time 
(t£ ) with respect to the barrier separation d for a wave packet having 
a narrow wave number distribution around the central wave number 
ko- All quantities are expressed in atomic units: 7 = 16, ko = 1.2, 
a = 20, zo = — 8cr. The mean dwell time t/j, the average reflection 
time (t™ V and the free time tf Iee are also shown for comparison. Fig 
(a) corresponds to the range of small to moderate barrier widths (in 
which we used a log scale for the vertical axis). Fig. (b) shows a wide 
range of barrier widths [we have multiplied all the times, except (t^), 
by a convenient factor (= 130), in order to display all the curves in 
the same plot]. 
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for the SWP clock. That and (if) are good candi- 
dates to transmission and reflection times also follows from 
the fact that they naturally satisfy the generalization of 
the "weighted average rule" [57] for the time dependent 
case [eq. (18)], provided the asymptotic conditions are 
satisfied. Also, the above approach has the advantage of 
providing these average values as the only time scale, in 
contradistinction with previous treatments of the Larmor 
clock in which they appear as the real part of a complex 
time [HI EH [28] . 

We analyzed the properties of the average transmission 
and reflection times for the scattering of a gaussian wave 
packet, centered on a tunneling wave number, by a square 
barrier and a double delta barrier and found that they 
would deviate from the mean dwell time only when the 
transmission or the reflection coefficients were very small 
[for symmetric potentials tj(k) and if (fc) equal the sta- 
tionary dwell time and the average times fl4[) and (15) 



are an average dwell time weighted by the transmission 
and reflection probability density, respectively]. Then, by 
carefully considering the opaque limit we found that (t£) 
shows no evidence of either the Hartman effect [7J (single 
barrier) or its generalized version [5] (double delta barrier) , 
and it also shows that the tunneling particle is slowed down 
by extremely opaque barriers. These properties indicate 
that (t£) is a good candidate for an average transmission 
time. It should also be noticed that our conclusions con- 
cerning the emergence of the generalized Hartman effect 
are in disagreement with those of HI] (who used a dif- 
ferent definition for the average transmission time) in the 
region of small barrier spacing, for they claim the exis- 
tence of such effect for a wave packet as long as the width 
Ak of the wave packet is much smaller than the spacing 
between resonances, while our results for (t^) show that 
the generalized Hartman effect does not emerge in any re- 
gion of the parameters. We stress that average tunneling 
times are the most one can expect due to the lack of a well 
defined time operator and the consequent impossibility to 
obtain an eigenvalue for time. 

The above analysis shows that the saturation of the 
transmission time, a feature associated with the Hartman 
effect, is a consequence of not taking into account the fact 
that the definition of a meaningful traversal time asso- 
ciated to a tunneling particle requires localizability, a re- 
quirement impossible to be satisfied in the stationary case. 
In addition, it is also a consequence of not taking the prob- 
ability of transmission into account properly. 

It should be noticed that even though the average trans- 
mission time does not saturate, for intermediate barrier 
widths (t^) increases very slowly with a (single barrier) 
or d (double delta), so that the transmitted wave packet 
"speeds up" in comparison to the free wave packet (see 
Figs. [I]and|4]). Although this phenomenon is well-known 
[4] , it could be interpreted as a "remnant" of the Hartman 
effect and it deserves further investigation (in the realm 
of relativistic quantum mechanics) due to the fact that it 
could allow for (non-causal) superluminal average speeds 



in this range of barrier widths. 

Finally, given the sensible properties exhibited by the 
average times {tf) and (t£), it would be interesting to 
consider its possible application to the analysis of quantum 
tunneling in graphene, a field of much current interest and 
in which many of the issues discussed in this paper are of 
practical relevance (see, e.g., [25] and references therein), 
including the controversy around the Hartman effect [HU 

ESQ. 
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